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ABSTRACT 

The  approximation  of  solutions  of  elliptic  problems  of  order  2 1  over  two-dimensional 
polygonal  domains  by  the  p-version  of  the  finite  element  method  is  investigated.  Op¬ 
timal  rates  of  convergence  are  established  for  the  case  when  elements  possessing  C<_1 
continuity  are  used. 
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t.  Introduction 


Thr  finite  element  method  has  three  versions:  the  h-version ,  the  p-version  and  the 
h-p  version.  In  the  /i-version,  increased  accuracy  is  achieved  by  decreasing  the  mesh 
size  /)  while  keeping  p,  the  degree  of  elements  used  fixed  (usually  p  —  1,2,3).  In  the 
p- version,  a  fixed  mesh  is  used  while  the  degrees  p  of  elements  are  either  uniformly  or 
selectively  increased  to  achieve  accuracy.  The  h-p  version  is  a  combination  of  both. 

The  standard  //-version  has  been  thoroughly  investigated  and  many  commercial 
and  research  programs  axe^ayailable.  In  particular,  the  solution  of  elliptic  problems 
of  order  2 (  using  elements  that  are  C*  1  continuous  has  been  analyzed  and  optimal 
convergence  rates  established  (see,  for  eg.  [8]  and  the  references  therein). 

:  The  />-  and  h-p  versions  are  recent  developments.  There  is  only  one  commercial 
code,  the  system  PROBE  TVHrpSiT'-IySt>ii)  "and  the  first  papers  discussing 

theoretical  aspects  appeared  only  in  198 li  ([6|,[2]).  For  an  account  of  todajPi^state  of 
the  art,  see  |l|.  See  also  |3|,  |4|,  [5|,  |9|,  jlO],  ( 12)  [13). 

In  [r>|,  a  second-order  mod'd  elliptic  problem  was  considered  over  a  bounded,  two- 
dimensional  polygonal  domain.  It  was  shown  that  if  (7°  elements  belonging  to  //„ 
were  used,  then  the  rate  of  convergence  was  optimal  up  to  an  arbitrarily  small  r  >  0, 

MU'  <  <W*  "'>11,1.  (i.i) 

The  dependence  upon  (  was  removed  in  (3)  to  yield  an  optimal  convergence  rate  for 
( (dements,  namely 

Mu'  <  Cp~[k  1,||tt||ff.  (1.2) 

Several  higher  order  problems  occur  in  engineering  however,  like  problems  of  plates 
and  shells  (|I9],  [20] )  for  which  elements  possessing  a  greater  amount  of  continuity 
than  elements  are  required.  To  this  end,  the  results  of  [6]  were  extended  in  [16] 
lo  ( ' 1  elements  to  obtain  the  following  rate  of  convergence: 

Mill,  s  (1.3) 

Unfortunately,  for  the  case  where  the  solution  lies  in  Hk  n  //„,  the  proof  in  [16]  is 
predicated  on  an  interpolation  assumption  ([16],  eq  (2.37))  without  which  (1.3)  does 
not  hold.  This  assumption  is  stated  by  us  in  Section  4,  Remark  4.2.  Moreover,  the 
proof  of  |I6|  indicates  that  the  term  (7(r)  can  grow  quickly  with  c  -  ►  0.  Nevertheless, 
computational  experience  indicates  that  (1.3)  holds  without  the  term  f. 

In  |9j,  an  alternate  approach  is  used  to  show  that  for  Cf  1  elements,  £  >  1, 

Ml „.<(.'(< )l>  “  '"•Il«lll,.  (1-4) 

The  proof  of  (1.4)  does  not  specifically  deal  with  the  case  of  boundary  conditions. 
In  order  to  extend  (1.4)  to  conforming  ('/  x  elements  with  f  1  vanishing  normal 
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derivatives  on  the  boundary,  an  assumption  similar  to  the  one  in  [16]  is  used  implicitly 
in  the  proof  of  [9],  Theorem  3.4. 

In  this  paper,  we  investigate  the  approximation  of  functions  in  Hk  n  //^  by  the 
/>- version,  using  polynomials  in  C'  1  fl  //*.  For  k.  >  2 t.  we  obtain  the  optimal 
approximation  result 

||ell/r  ■<J||u||w»  (1.5) 

In  proving  (1.5),  we  do  not  use  the  interpolation  assumption  used  in  |9],  [16].  The 
use  of  this  assumption,  however,  allows  us  to  extend  (1.5)  to  the  case  (  <  k  <  2f 
In  the  case  for  which  the  solution  exhibits  singular  behaviour  of  the  type  v  rn 
((r.O)  being  polar  coordinates)  and  the  vertex  of  the  elements  is  at  the  origin,  we 
j  obtain  the  optimal  estimate 

I 

IMI/r'  <Cp  (1.6) 

This  improves  upon  the  rate  of  convergence  found  in  [16]  (and  |9]-[l0|)  which  is 
optimal  up  to  an  arbitrarily  small  t  >  0. 

In  Section  2,  we  describe  the  notation  used  and  our  model  problem.  Section  3 
deals  with  approximation  properties  of  polynomials  on  the  square.  In  Section  4,  we 
analyze  the  rate  of  approximation  for  functions  in  n  //„  and  prove  (1.5).  Section 
5  deals  with  the  case  when  the  solution  exhibits  singular  behaviour  and  proves  (1.6). 
In  Section  6,  we  summarize  our  results  and  briefly  address  some  generalizations. 


2.  Preliminaries 

2.1.  Notation 

bet  R2  hr  t ho  usual  Euclidean  space  with  x  —  (xi,.r2)  f  R2.  fl  C  R2  will  denote  a 
hounded  polygonal  domain  with  vertices  /!,,  t  ~  0, 1 , . . . ,  M,  A0  ~  AM  and  with 
boundary  I'  —  T\  where  P,  are  open  straight  lines  with  end  points  The 

inlernai  angle  between  P,  and  PH  (  will  be  denoted  by  w,,  t  =  t, . . . ,  M,  0  <  w,  < 
2 7i.  'Phe  case  u,  —  2n  results  in  a  slit  domain  for  which  the  boundary  is  two-sided  (in 
an  obvious  sense). 

My  7/2 (n)  —  fV\n)  and  flk( H),k  >  0  we  denote  the  standard  Sobolev  spaces 
(with  index  2).  Also,  //* (11)  denotes  the  subspace  of  functions  with  k  -  1  vanishing 
normal  derivatives  on  P.  For  k.  >  0  not  an  integer,  we  define  //*(fl),  //*(fl)  as  the 
usual  interpolation  spaces  (by  the  A-method,  see  [7]): 

//'*'(  n)  {fi'(n),  n'"  {n))fi<l 

where  f  \  0  k,  0  <  0  <  1 ,  q  -  2.  For  k  >  1 ,  we  define  /f*4*(n)  =  //*+*( n)  D 
llk(U).  We  will  also  deal  with  the  Sobolev  spaces  /7*(P,),  77* ( /) ,  /  =  (n,  b)  which 

are  defined  for  k  integer  in  an  analogous  way.  The  spaces  //*(!!),  //* (n) ,  //*(P,), 

etc.  are  Hilbert,  spares  and  their  inner  products  will  be  denoted  by  (•,  *)//*(n)i  etc. 

For  k  ■  0,  we  let 

tt(K)  {(T|..r2)  |  |x,|  K,  1.7*2 1  <  (2.1) 

///%;«( ll(*))  1  77* ( /?(«■))  will  denote  the  spare  of  all  periodic  functions  with  period 

2k. 

Mv  P,!(ll),  respectively  7J,1  (/?(«■)),  we  denote  the  space  of  all  algebraic,  respec¬ 
tively  trigonometric  (with  period  2k),  polynomials  of  total  degree  at  most  p  on  11, 
respectively  R[k).  My  P*(ll),  respectively  T2(R{k)),  we  denote  the  spare  of  all  alge¬ 
braic.  respectively  trigonometric  polynomials  of  degree  at,  most,  p  in  each  variable  on 
If.  respectively  R(k).  We  also  define  fP(P,),  Pr(l)  (/  —  («,/>))  to  be  the  space  of 
all  algebraic  polynomials  of  degree  at.  most  p  on  P,,  respectively  I. 

bet  1?  i  ,11,  where  11,  are  (open)  triangles  or  parallelograms.  We  shall  assume 
that  {],  i  i  11,  -  <f>  for  i  /  j  and  11,  n  11;  is  either  the  empty  set  or  an  entire  side  or  a 
vertex  common  to  11,  and  11;.  We  will  assume  that  all  vertices  of  11  are  the  vertices 
of  some  Uj.  I!,  will  he  called  elements. 

I<et  Q  (  1,1)  (  1,1)  anrl  T  {(.r,,T2)  ||.T||  <  1,  I  <  ,r2  <  r,}  denote  the 
standard  srpiare  and  standard  triangle  respectively,  bet  l\  be  an  affine  mapping  with 
.buohian  having  positive  determinant  which  maps  II,  onto  Q  if  11,  is  a  parallelogram 
and  onto  T  iT  11,  is  a  triangle,  bet  now  .^(U)  <  l2(U)  denote  the  set  of  functions 
»/  such  that  if  »/fl  denotes  the  restriction  of  u  to  11,  then  Un ,o(l\)  1  r  P2{Q)  if  11, 
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is  a  parallelogram  and  wno(/'’,)  1  f  Pp{T)  if  H,  is  a  triangle.  We  will  then  write 
Mu,  '  *>(11,)  an(l  ‘  M11)-  Furthermore,  we  denote  for  t.  >  1  integer, 

•‘v'W  -  //'(fl)n  sp([]) 


oS'(u)  s'(n)  n  H'( n) 


It  is  possible  to  show  that  S'  >  ( '!*  l>(f2)  and  uSp  <  (%  ’’(H),  where  C^(Ct)  is  the 

spare  of  all  functions  with  f  rontinuous  derivatives  on  fl  and  C (_  C^(Cl)  the 
suhspaee  of  functions  with  compart,  support  on  fh 

For  r  ( r ! ,  r2 ) .  ()  r,,r2  |r|  r,  I  r2  and  u  a  function  defined  on  fi,  D^u 

will  denote  the  derivative 


l){r)it 


u 

f?T|'  c?t23 


2.2.  The  Model  Problem  and  its  Properties 


l«et  /,  be  an  operator  defined  on  '  I  by 


Lit  (  l)'A'i»  I  (  I) 


'  Ia'  I 


u  . . .  An  I  u  -  y  (  I  y  A1  u 

,-a 


(2.2) 


I, el 


It} 


denote  differentiation  with  respect  to  the  outward  normal  to  fl  on  F.  Then 


wo  will  consider  the  following  model  problem 


hn  —  F  on  n  (2.3) 

iYu 

q  "0  on  F  r  0, 1  ~  1  (2.4) 

Note  that  (2.4)  is  equivalent  to 

l)(r]n  -  0  on  F  0  <  [r |  >  f  -  J 

The  model  problem  (2.2)-(2.4)  is  a  classical  case  of  an  elliptic,  equation  over  a  non- 
smooth  domain.  The  properties  of  such  problems  have  been  well  studied  (see  (17), 
||N|  and  the  references  therein). 

We  assume  here  that  the  solution  of  (2.2)-(2.4)  can  be  written  in  the  following 

form: 


u 


M 


n  i 


“2 


(2.5) 


5 


whore 


u,  e  //*(n)  n  //'( n) 


4'1  -  RejXjcf  llogr.l^"^  |  r 


jo 


(2.6) 


with  R.e(o }  )  '  f.  1,  Re(n^'] ,)  >  Re(o^),  7^  >  0  real,  4>^{0t)  and  x^(rf)  are  real 

t"'  (or  sufficiently  smooth)  functions,  x^(ri)  ~  1  for  0  <  r,  <  pl'l  <  ,  x^(r<)  —  0 

for  r,  >  2 pi'l.  By  (r, ,  0,)  we  have  denoted  the  polar  coordinates  with  the  origin  at 
the  vertex  /I,  of  the  polygon  fi.  The  partition  (2.5)  is  typical  for  problem  (2.2)-(2.4), 
with  the  functions  t/.j^  describing  the  singular  behaviour  of  the  solution  caused  by  the 
corners  of  R.  For  details  and  proofs  of  the  partition,  we  refer  to  [17],  [18]. 

It  may  be  noted  that  we  have  only  dealt  here  with  essential  homogeneous  boundary 
conditions.  Instead  of  (2.4),  we  ma.y  specify  natural  boundary  conditions  as  well 
(which  may  be  inhomogeneous),  consisting  of  normal  derivatives  of  order  i  <  r  < 
2(  1.  Our  results  remain  valid  for  the  case  when  different  types  of  conditions  are 

specified  on  different  portions  of  I'.  In  section  6  we  shall  comment  on  this  and  other 
generalizations  of  our  results. 

2.3.  The  p- vers  ion  of  the  Finite  Element.  Method 

Problem  (2.2)-(2.4)  may  be  put  into  the  following  equivalent  variational  form:  Find 
n  f  //*(R)  satisfying 


/?(«.,  t/)  -  /  FvdU  for  all  III 

J  o  J 


(R) 


where  f?(-,  •).  the  bilinear  form  defined  on  //'( R)  x  Hf{U)  associated  with  the  operator 
//  is  equal  to  (vjff'fn)  so  that 


(2.7) 


»(«,«)>  IHlir'tfi)  (2-8) 

holds  for  any  u  £  Hf(U). 

I’he  p-version  of  the  finite  element  method  consists  now  of  finding  ur  £  0^(0) 
such  that 


n(nr,v) --  /  FvdU  forallwf  0^(H) 

J  r»  " 


(2-9) 


B.v  the  results  of  [8|,  for  p  --  4f  3,  n <‘>p  consists  only  of  the  function  0.  Hence, 
using  (2.8),  (2.9)  gives  us  a  unique  riv  f  0  for  all  p  >  4 1  3. 


mmm 
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3.  Approximation  Properties  of  Pr(Q) 

I, rl.  Q  hr  (hr  standard  square  with  sides  7,,  i  -  1 , 2, 3, 4  and  diagonal  75  corre¬ 
sponding  to  Z|  -  r2.  Let  T  he  the  standard  triangle  enclosed  by  sides  7, ,72  and 

7v 

Lemma  3.1  Let  S  Q  or  T.  Then  there  exists  a  family  of  operators  { ft,, } ,  p  = 
1,2,...,  IIr  :  Ifk{S)  ->  Pr(S),  such  that  for  any  u  6  Hk(S), 
for  k  >  0,  0  <  <7  <  k: 

||M  IIPMII //»(.<?)  5-  ('P  **  7^IIM II //*■(.'?)  (3-1) 

for  r  -  ( r , ,  r  2 ) ,  k.  >  |r|  t  0  v  s  <.  k  |r|  i  —  1,2,...,  5, 

\\l){r](u  Hr«)ll/f*b.)  <  Cp ■<* >' *-  i»||ti||,f.(5)  (3.2) 

/or  r  -  (r|,r2),  k.  >  |r|  It,  re-  .S', 

|P(r)(«  Urw)(.r)|  <  C/»-(t'|r|-,,||tt||H»(5,  (3.3) 

where  the  constant  C  is  independent  of  11  and  p  and  where  we  denote  /)r(><’)  ~  Pp{Q) 
for  S  Q  and  FP[S)  -  P}{T)  for  S  -  T. 

Moreover,  if  ti  C-  /^.(iS),  then  IIr»  u. 

Proof  -.  Let  d  >  \.  Then  S  c  R(d)  (see  (2.1)).  Since  S  is  a  Lipschitz  domain  there 
exists  an  extension  operator  T  mapping  Hk{S)  into  Hk[R{2d))  such  that. 


T11  -  0  on  R(2d)  -  R{\d)  (3.4) 

where  <7  is  independent  of  u.  A  concrete  construction  of  T  may  be  found  in  [21). 

Let  <h  be  the  one-to-one  mapping  of  R(2)  onto  R[2d): 

R(2d)  ->  x  -  (j-,,.t2)  -  -  (2dsin  <£,,2dsin  £2) 

with  (£,,£2)  £r  H{\)-  hurt  her,  we  let 

R  '\R(ld)\c  R(’) 

where  <f>  1  denotes  the  inverse  mapping  of  <f>. 

Let.  v  T11  11ml 

VU)  -  »W))-  (3-n) 


‘jQSk. 


Because  of  (3.5),  we  easily  see  that 


supp  V  (£)  C  R  (3.7) 

Tn  addition,  it  can  be  readily  seen  that 

(a)  V7(£)  is  a  periodic  function  with  period  2n 
(l>)  V  satisfies 


li*7  ( 0 II »*( id*))  -  CIM!  nvW\d))  S  Cllwllz/Hs)  (3-8) 

and  hence  V  6  ^e/?/?(^(7r))- 

(c)  F(£)  is  a  symmetric  function  with  respect  to  the  lines  £,  =  ±y,  t  =  1,2. 
bet,  us  expand  the  function  V  in  terms  of  its  Fourier  Series 

CO  OO 

b'(6,6)  ‘  E  E  (3.9) 

J  •  OO  f.~  -  OO 

For  any  />  ;>  I,  we  define  Vr  —  flp  V  by 

M£i,6)  -EE  «*e<(,'€,+*a)  for  S  —  Q  (3.10a) 

i~-v  t--r 

^(6,6)-  E  «i<e,'(i€"7t4)  for  S  =  T.  (3.10b) 

l;IH'l<r 

Obviously,  V7r  Tr(R(n)). 

We  have 

■>  ? '  (3.U) 

j.t 

where  has  the  usual  meaning  of  equivalency.  (3.11)  yields  immediately  for  0  < q  <  k. 

II’7  KWnum*))  '  cp 

<:P  (*  ,)ll»llnno  (3.12) 

using  (3.8b)  with  ('  independent  of  u. 

In  what  follows,  we  assume  S  Q.  The  rase  S  T  follows  similarly. 

l-et  7,.  >  ~  1 . 4  be  the  sides  of  R(i r)  and  let  £2  £.2  be  one  of  the  sides.  Then 

fm  r  (r, .  r2),  0  jr|  -  k  J, 


4 

4 


1 


D^(V  -  Kr)(6,6) 


\|j|>rM<p  bl<rl'l>p  !;!>;>  \(\>r) 

-  21  i  E(or4V"'  f 

l>l>r  |>|<p  |>|>p 


where,  for  |/|  >  p: 


"Y'h'P |2  -  fr,{}2  tr7ajeeie^)2 
M<p 

i  !’>!r’(E  M)! 

I4<r 

.?'2Vr3( E  M20  +  J2  +  /*)*)(  E  0  +>2  +  «2)-*) 

ki<p  bi<p 

C  j2r>  p2r*  * 1 A  j 

(11  f)r 


where  we  denote 


E  M*0  if  +  '*)*• 


Consider  now  t. ho  fund, ion 


for  r  f).  VVp  have  for  //  -  0 


/(*)  - 


(1  I  .T2)* 


(3.H) 


f'(X)  '  (1  ,  r2)*  *  I 

IIpiu  p,  for  k  >  //  ami  :r2  >  j/‘  ,  /(j)  is  a  decreasing  function  of  ,r.  Moreover,  when 
I1  H,  /(.r)  is  decreasing  for  all  .r  -  0.  Hence,  (.here  exists  a  constant.  C  —  C{k,fi) 
such  that  for  |.;|  -  p,  k  .■>  //, 


72"  Cp2" 

(I  t  72)*  ’  (1  I  P2) * 

so  t  hat  I  aking  //  r* ,  for  A-  r! , 

k«7T,H,,i*  •  ,2lr  2|r|  ” 


(3.1  -V 


(3.16) 


**Vi”*v 


K’jY'b  12||2 


1*1  >r 

**■>»>(£  ^0  +  /  +  «2r‘) 

i*i>p 

dx 

■^;>  +  l  X^(k~rA 

CAjP~{  n-2"-" 


Cp2r'  A}  /' 


provided  k.  >  r2  f  J. 
Analogously,  for  |j|  > 


|(,:,)W,fV  <  ,*• 


CM 


Vl*l» 

oo 


,1 


j2r,X2r‘1 


<  CM 


p+i  (1  +  j 2  4-  x2)* 

1 00  dx 


dx 


pH  (1  +  j2  +  XZ)<!'r,~r2 

CMj7r<2*  2H-» 


provided  /r  >  j r | 


lienee,  using  (3. 13)-(3. 18),  we  see  that  for  i  =  1,2,1 


ll^,r)(v  -Wl/mj  <  c 

<  cP £  yj. 


£  l/'fty  l2+  £  l’+  £  I/'6! 

I>l>p  lil<r  l;l>r 


;“-oo 


<  CM  2<t  |rl  ■sMl.ill* 


P  ‘  ll"llff*(Q) 

provided  k  >  jr|  f  ‘ .  From  t  his,  it,  follows  easily  that  for  s  >  0,  ,s  integer, 


ll^(r,('/  vp)||,/.m<c>^-w-‘-ii||»t||/f»W) 

provided  k  >  ] r |  |  .s  I  ^ . 

Now  we  estimate  ||f4(r)(l/  ^lOllwpr-)-  For  ^  ^  |rl  ^  k  ' 


'  I  .1  ^  I  i'l  _  I  .1  _  I  tf  I -  I.  I- _ I  c\ _ ' 


IjI  >»■  Kl<r  IjI<p|£|>p  hhr>M>r 

OO 

X  (<•'.!''  I  f.f  I 


e!"  -  E  (<»'■  (>«)'■  <* 


J  I  '~7 

hl>r,k|<r 


'f  -  E 

;  t  '-7 
URr.l'br 


Cf  -  E  ('»ri  (><)”«,/ 

f“7 

|j|>p.|P|>r 


Ry  Mi r*  Schwarz  inequality, 


j  t 

,1>p.K!<r 


E  M2(<  <  j2  +  m[  E  >2r,^2rj(i+y2  +  f2y 


;  <  t~q 
hl>r.|*l<r 


vj" 


■vhere 


/U/  -  E  lvl2(<  I  .?2  1  'Y 

;  I  f  -7 

/i’1  -  E  (<7  ■  02r,^r3(l  I  .7 2  i  /"2)“* 

j  t  (-? 

IjI  r-|p|<f* 

'  22r'  1  E  (f/2r'  1  'c2r,)^2rj(1  I  J2  I  **)  * 

;  i  f  7 

b'br.l'Rr 


''or  q2  it 2 ,  this  gives 


r 


where  !V(p,q),  the  number  of  terms  in  the  sum,  is  obviously  <  2 p,  so  that 

/lb  <  c>  (2fc -2|r|-l) 

For  p2  72,  we  have 

<1 2  -  (./  i  q2  <  2 (7 2  1  t2) 

so  that 

<  Cp  (2*-2M  ’) 

provider]  k  >  |r|.  Hence, 

|C’],||2  <  CAqp  {n  2H  1) 

Similarly, 

|Cf|2  <  CM,p-'2^2H-» 

Finally, 

|cf  |!  <  r,M,/i’i 

where 

; 1 2 1  v-  (q2r'  \  P2r')f 

7  "  }7tq  (1  kj2M2Y 

\A>rM\>r 

i  •  2  2 

ror  q  '  P 

,|2|  y  (?’2r’  l  P21')^ 

ik  (i  •  nk 
<  r-n  (2*  2|r|  1) 

provided  k  >  jr|  f  ^ . 
If  q2  "  p2,  then 


/I2) 

7 


r  M-'j 
bl>r.l'l>r 


q2r'(2r' 

(i  »  p  t  f2y 


v  '<  . _  ,  pn-iH-jirM) 


12kl 


J  *  >'7 

Id  'r.r"  M-'M 


S  C’|7 
<  Cp'(w  ■2H'1) 


q‘<'  v  ,  f2lr 

2w 

l'l>kl 

(**•  2H-1)  +  (7p-(2‘  2M  ') 


V  V  I-  \ "  <  r/n  (2*  2H- 1) 

^  (i  +  v>  Mru,(*»^)‘ 


provided  A:  >  |r|  f  ^  Hence 


|6f  l|2  <  CAqp  <2t  2H  » 

Combining  (3.20)-(3.23),  we  obtain 


(3.23) 


l|T>(r)(T  -  ^ll/Hi,)  <c>"(2*-2lr|-,) 

which  leads  to  (3.19)  for  7r,. 

We  now  estimate  | D^(V  -Cp)(x)|.  Recause  V  -  Vp  f  HpFn(R(ir)),  we  can  assume 
without  loss  of  generality  that  (6,  £2)  C  7,  given  by  £2  ~  Ci-  Then  once  again  we  have 
(3.13)  and  for  c  >  0 


Taking  // 


r  1  I  |  -I  (  in  (3.l.r>),  we  obtain  analogously  to  (3. 10) 


^CA,p 

provided  k  >  r,  -(-  |  I  f,  so  that  (3.24a)  yields 

\bl»>  / 

Moreover,  using  (3.17),  we  see  that 


(3.24a) 


(3.24b) 


ijd. 


d;Kr 


\bl<r  / 

r1  2(<t  . 


13 


provided  k  >  r2  I  ’.  Also, 


'EKor^l)  <  (e  -t!*) 

Vbl>r  /  llj|>r  /  \W>r  •*  ) 

<  c  f  Y.  AiP  {2k  2|r|  2~2f])  p  *' 

Vw>r  ) 


<  Cp -a(*-|^,)||H||J#M9) 


(3.26) 


provided  k  >  jr|  +  1  I  t,  where  we  have  used  (3.18)  with  r\  =  rj  +  =  f  £.  Combining 
(3.21b)-(3.26),  we  get  for  k  >  |r|  f  1, 

\f>lr)(V  Vr){x)\<CP  ,fc-|r|  ,>||»/.||/,ng)  (3.27) 

Because  of  (3.8r).  Vp(<f)  '(t))  f-  Pr(Q).  Further,  <f>  is  a  regular  mapping  of  R(d.) 
on  Q,  (d  r  * )  and  4>(li)  >  h-  Hence,  for  k,q,s  integers,  (3.1)  follows  immediately 

from  (3.12),  (3.2)  from  (3.19)  and  (3.23)  and  (3.3)  from  (3.27).  Using  a  standard 
interpolation  argument,  (3.l)-(3.3)  will  hold  for  non-integral  values  of  k,q  and  s  as 
well. 


The  following  one-dimensional  result  is  from  |5|. 

Lonmin  3.2  I,et  I  (  1,1)  and  u  (7.  k  >  1.  Then  there  exists  a  polynomial 

zv  ‘  p  >  2k  I,  such,  that 


z},'}(  I  I)  -  u(,)(  I  1)  t  -  0,  1 . k  -  I 

and  for  s  —  0,  1 , 2 . k. 

Il»  *rllw'(f)  <  CP  [k  f) II” II //*(') 

where  ( '  depends  on  k.  but  not  on  u  and  p. 


(3.28) 


(3.29) 
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4.  The  Approximation  of  Functions  in  //fc(l?) 

In  this  section,  we  will  analyze  the  rate  of  approximation  in  the  //'(fl)  norm  of  u 
by  a  piecewise  polynomial  in  0 5^(11)  for  the  case  when  «  6  //*(!!)  n  //^ (fl)  ie  when 
wl'1  0  in  (2.5).  The  main  result  we  obtain  is  the  following. 

Theorem  4.1  Let  f  be  an  integer,  (  >  l.  Let  u  6  //k(fl)  n  k  >  21  —  j. 

Then  there  exists  a  sequence  { r/. r } ,  p  >  4f  -  3  satisfying  up  G  oSp(ft)  onrf 

^0  r  —  0,  I  -  1  on  F  (4.1) 

II"  Mrll//'(n)  ^  C/*  ,V  /*llull//‘-(n)  (4-2) 

where  the  constant  C  depends  upon  the  partition,  of  fl,  k  and  f.  but  is  independent  of 
ii  and  p. 


We  first  prove  the  following  lemmas. 


Lonmin  4.1  Let  /  —  (  1,1).  Given,  an  integer  1  >  1,  there  exists  a  sequence  of 
polynnmi  *  {Xp}  ^  {.Yp'h  P  >  2f  I,  r  0,1,...,/  ],  in  Pr(J)  satisfying 


dr1 


(  I)  -  I  «/  7  -  r  j  -  0,  !,...<  1;  r  -  0,  1 , . . . ,  t  1, 


^x; 

dr 1 


n  if  j  f  r 

II)  -  0  -  0,1 _ f  I;  r  -  0,1 . /  -  1, 


(4.3a) 

(4.3b) 


ll\pl|//*(/)  ('P'  '  ’  *  0,1 . t  I;  r  -0, I, 

where  the  constant  ( depends  on  t  but  is  independent  of  p,  r  and  s. 

Proof :  We  first  define,  for  p  '  2 1  l,  functions  <f>rv  given  by 

*;.(t)  I  l)’)r 

'  I  n  / 

where  the  constants  ( ’  are  given  by 

C,  r  0  for  0  <-  i  '  r  1 

l (!'' f°r  r  1  1 


(4.4) 


(4.5) 


(4.6) 
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I  loner. 


ICrl  <  Kpi  r 


By  (1.5),  for  any  integer  0  <  m  <  t  1, 


(4.7) 


(jin  ±T 


i-n 

<  l 


}  m(-p)m  1)‘  V  *'♦»  f 


i- | 


f  V  r:r'p,!(x  1  1)'  "v  p(*,l) 
'  (*'  -  »«}! 


(4.8) 


so  that, 


dm6r  m  1  <  /  f  i  \  > 

h'p.flhnit  <  «E;>"  'l>ri  (I  (r  I  ip  »« 

>-n  t-j  '  >  / 

Now  for  p  sufficiently  large, 

/’V  I-  iYe  ^^dx  <  *><'+*) 

Hence,  using  (4.7), 


dm (f>r  ™  «-  i 

II  (lTJ IrrV)  < 

;  -o  i~j 
<  Kpm  r-  5 

libs  shows  that  for  any  s  —  0, 1 , . . . , <  -  I, 

ll^pll//-(/)  <  Cp"  r  i 

ir.  {,/.'}  satisfy  (4.4)  with  replaced  by  <fyrp.  We  now  show  that  with  our  choice  of 
r'.r',’>  a,so  satisfy  (4.5a).  By  (4.8),  we  have  for  0  <  m  <  t  -  1 


ArJ 


A  j  n  ./• 


I 'sing  (4.5),  we  get 


ArJ’  0  for  m  '  r  I 
I  for  m  r. 


Also,  for  r  <.  m  <  t.  1,  we  have 


■"  [ZUn-r-Mj 

since  the  term  inside  the  brackets  is  the  binomial  expansion  of  (1  1)m‘r. 

In  order  to  obtain  a  function  satisfying  (4.3b)  as  well,  we  let  U'  —  <f>Td)  where 
xj'  <-  (*'Y'(/)  is  a  smooth  cut-off  function  satisfying 


V’(x)  =  1  for  1  <  x  < 


0  for  -  <  x  <  1 . 

2  ~  “ 


Then  it  may  be  easily  verified  that  Iff  satisfies  (4.3)-(4.4).  We  now  use  Lemma  3.2 
to  approximate  Iff  by  a  polynomial  zrp  in  F,  (/),  p  >  21  —  1,  satisfying  (3.28)-(3.29). 
(3.28)  implies  then  that  zrp  satisfies  (4.3).  Also,  using  (3.29)  and  the  fact  that  Uf 
satisfies  (4.4),  we  have 


\\v;--*;\\ !/-(/)  <  cp-(fc-')iit/;iiw*(/) 


so  that  by  the  triangle  inequality  ,  zf  satisfies  (4.4).  Taking  \Tr  =  yields  the  lemma. 


Lemma  4.2  Let  v(x)  be  a  function  defined  on  I  satisfying 

d’  v 

.-.-;( I)  -0  s  ~  0,  1 , . . . ,  r  1 


(x  I) 


C lie  HO 


(4.10) 


I  Woof :  Let  t(£)  be  a  function  defined  on  [0,oo)  such  that  /  0  for  x  '  2.  Define 


/,(£)  '{fit  <)''/(»)* 


(4.11) 


'I’hen  bv  |I5|,  page  245,  No.  (9.9.5),  we  have 
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Now,  take 

/({)  "'"(I  £)• 

Then,  substituting  in  (4.1 1),  integrating  by  parts  r  times  and  using  (4.9),  we  see  that 

fr(  0  =  <?”(!  0 

so  that  (4.12)  becomes 

Substituting  1  —  £  —  r  gives  (4.10). 

□ 


I ’roof  of  Theorem  j.l 

l-et  II,,  ;  1,2.  .  .  ,  A'  be  the  elements  of  the  partition  of  12.  We  first  construct 

functions  ?j,'j  11,,*/^  as  in  Lemma  3.1.  The  lemma  is  applicable  because  a 

linear  transformation  maps  the  parallelogram  or  triangular  element  onto  the  standard 
square  or  triangle,  preserving  the  polynomials.  Hence  for  f  ->  1  integer,  k  >  f. 


It1'  2jo lln'tn.)  '  -  ('P  *,||«IIhh».)  (4.13) 

Our  first  step  is  to  add  a  function  f  Pr(  12,).  p  >  4 1  3  to  so  that  the  function 

?!'1  1  ?/],'*  satisfies  (3.1),  (3.2)  and 

/2(r,4'1  l){r)u  for  0  <  .  |r|  <  2t  2  (4.14) 

at  the  vertices  of  fl,.  bet  us  first  assume  that  f2,  is  a  parallelogram,  which  we  may 
take  to  he  the  standard  square  Q  without  loss  of  generality,  bet 


/2,r)(» 


I, 


whe  i  e  r  ( r  t .  rj ) .  Define  the  function 


»'l  L  n*r  1  \  p'  (  T>  )  (  r2  ) 

,r|  2/  2 


(4.13) 


(4.10) 


Here,  ,  i  1,2  are  as  in  bemma  1.1,  with  t  2 f  1,  so  that  i.  Pr(/), 

p  If  3.  This  implies  that  we  may  construct  i/q  >  PV(Q)  provided  p  ■  if  3.  We 
s rr  t  lien  t  hat  with  k  if  I , 
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t 


iiM,iiiff'(n.)  ^  E  E  in'(r)nix)r,'ii//'»(/)ikp,ii//'2(/) 

0<|r|<2<  2  tyUi-t 

<  E  E  p‘‘lW+V''  ,-»P<>-r,-*iHiw.(«w 

0<M<M-2  'i4*2<r 

<  Q>  “  H-17) 


whore  wo  have  used  (3.3)  with  k  >  |r|  f  1  io  k.  >  2f  -  1  and  (4. 4).  Moreover,  on  the 
side  I’ i  {(.r,  1)|  I  -  r  -  1},  we  have  hy  (3.3)  and  Lemma  4.1,  for  j  —  (j\  ,ji), 
s  •  0  integer,  0  (j (  1  s  <  2 (  2,  k  >  2(  1, 


E  !«(r,l  II 

(l<|r|<2r  2 

E  l"(r,J1 

0<r,(jj<2  /  2 


~dxf  dxf{~inH’(r' 


>1  II 


v; 

d.r{' 


iH'(Pi) 


<:p (t  n  ^  v*  n  * iiwii«*(fi,) 
CP  <*  W-'-HWnMn.) 


(4.18) 


(4.18)  will  also  he  true  for  the  side  r2  -  {(-  1 , y)  |  1  <  y  <  1 }  and  will  hold  trivially 

on  the  other  two  sides  of  Q ,  where  wt  —  0.  We  can  repeat  this  construction  for  each 
of  I  lie  four  nodes  of  Q  to  obtain  w} ,  j  —  1,2, 3, 4.  Then  defining  yj,'!  —  £>-i  ii’j,  we 
see  that  rj,'l  I  i/j,'!  satisfies  (3.1),  (3.2)  and  (4.14). 

For  11,  a  triangle,  we  assume  that  II,  is  the  rotated  standard  triangle  T  with 
vertices  /’,(  I,  1),  /2(1,  l)andP3(  1,1).  Define  o*rl  as  in  (4.15).  Let  p>  I  Of  6 
so  that  ft  |(p  2 f  I  l)/2|  '  4 f  3  and  define  <-  Pr{T)  by 


1  *,)»'(  -2)  (1.19) 

with  \|!  '.  Then  we  see  that  for  0  '  |r|  <  2 f  2,  k  >  1, 

D(r,»r,,n°)  (Ion  fy,, 
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We  now  define  u< 


b) 

i 


recursively  for  |r|  >  0.  Lot. 


(-) 

(m) 


4 


H 


will  l>o  non-zero  only  when  r, 


1,  1).  (4.20) 

m„  t  1,2.  Now  defining  ft>r  0  '  |j|  <  2f  2 


»>(,;)  ("0)  X!  W)(  2)  Ir,)1'  1  (4.21) 

HM>I 
(* - ; .2) 

(wit.li  \j'  \^f  '),  we  see  that.  1,  I)  —  0  whenever  r,  <  j,  (t  —  1,2)  and 

i;1(  1.  -  1)  ~  f"01  E  |m|<|>|  /*&’)  for  r  -  j.  Hence,  for  p  >  \0f  6, 

'  n<  (t-  1,2)  7 

MM  -  1Z  «'i>)  b  M»)  (4.22) 

satisfies,  for  0  <  |r|  |  ,s  <  21,  2,  A:  >  2t  -  1, 

/;,r)te,(  1,1)  ^  ry<r)  (4.23a) 


D[r)wx  0  on  P2P,  (4.23b) 


ll«M||/i'(n.)  b  Cp  (*  ^Hullz/nn.) 


(4.23c) 


||/l,r,»>i || //.((■)  <  Cp  <*  ■|r|-,-i,||tt||ff*(n.)  (4.2.3d) 

wliere  f  f\  f\  or  /  2/V  lror  4 (  3  '  p  '  I0f  6,  we  may  still  construct  a  w j 
satisfying  (4.23a.b)  using  the  results  in  |8],  Then,  by  suitably  adjusting  the  constant 
C  in  (1.23c,d)  it  is  seen  tht  »/>t  ‘  Pr(ll,)  satisfying  (4.23)  may  be  constructed  for  any 
/>  \f  3. 

A  similar  construction  as  that  of  wx  may  be  used  to  obtain  functions  i/>2  and  te* 
associated  with  the  nodes  f\  and  l\  respectively,  after  first,  mapping  11,  in  a  suitable 
way  onto  the  standard  triangle  T.  Then  we  take  yhl  -  ( 

bet  now  7  Ut  i  i  Um  and  At.  A2  be  the  end  points  of  7.  Now  fVr)~W  f  fib) -HI 
on  Denote  the  “jumps"  of  on  7  bv 


I  hen ,  because 
/  1.2:0  |  r  | 


M 

jn< 


(r.Vl 


for  (  r.  y)  ‘  7 


P<r,4;|(d.)  D{r)v(A,)  /f,r|4m|(.4,),  we  have  w-JrJ(/1,) 

If  2  and  also,  by  (3.2).  (1.18)  and  (1.23d), 


(1.21) 


0  for 


for  0  | r |  (  /  <  2f  1,  k  >  2 ( 

Let,  now  Fj  ho  the  affine  transformation  satisfying  F^flj)  —  S,  where  S  —  Q  — 
{(£,y)  |  |£|  <  1,  |y|  <  1 } ,  the  standard  square,  if  Q;  is  a  parallelogram  and  S  -  T  = 
{(s.f/)  I  |£|  <.  1,  ■!<»/<.-£}.  the  rotated  standard  triangle,  if  U}  is  a  triangle. 
Let  7  he  mapped  by  F}  onto  the  side  I'j  —  {(£,  1)  |  |£|  <  1}  of  S.  Let  n  —  n (£,??) 
denote  the  outward  unit  normal  to  5  along  dS  and  let  n  he  its  image  on  dll,  under 

/■'  1  sri  I  hat  ,?-f  I  -  ,5/ 

1  j  ■  so  ,  ha'  at,  I’ll,  an 
integer. 


|as  whenever  /(x,y)  -  /(£,  r/).  Define,  for  (r,y)  €  7,  a  >  0 


ft"  (r,  y)  -  (zM  -  2|mh  . 

V>  0n'[p  p  V 

'•’lien,  if  li’m  0'moF}  1  we  see  that  for  (x,y)  e  7,  (£,  1)  =  (^, y) , 


(4.26) 


(4.27) 


\,nv  ,y;m 
satisfy  |r| 
2f  I,  Hr 


are  obviously  linear  combinations  of  those  w  defined  hy  (4.24)  which 
s.  Hence,  it  can  he  deduced  from  (4.2.r>)  and  (4.27)  that  for  0  <  s  |  (  < 
'  2f  ! , 


ll^,mll/t'(r,)  "  ('P  1  *,dl«Hw*(nJ)  *"  llull//t(n,„))  (4.28) 

Moreover,  since  for  |r|  -  2 f  2,  /J^zj/I  -  at,  the  end  points  of  7,  we  have 

for  I  f\  I , . . . ,  2 f  2  .s. 


d'/f 


Lor  U}  a  parallelogram,  define 


(")  " 


t  1 


(4.29) 


s  jtn 


(£•>/)  IXjovM 

n 


wifh  \*  \*  .  Thru  using  (1..X),  (1.20)  and  (f.2*),  vvr  s rr  that 

{j"‘  iF  on  L,  for  O  s  f  I. 

oil  ’ 


(4.30a) 


D,rl  C„ 


U  on  itS  I’|  for  ()  '  |r|  •  f  t 


(1.30b) 


f  1 

IlCjm  ||//'(.C)  ('  Yi  H^>»n  H/l/l(r,)ll\fll//'l(/) 

(,  \t?<  t  r-n 

<  c  Y  X>  “  '  i||«||ff*(n, 

/,  t  t2<t  .«~0 

<  <*  '>||u||f#*(n)  (4.30c) 

where  we  have  taken  A-  -  2f  I. 

We  now  show  (lie  existence  of  a  <jm  satisfying  (4.30)  in  the  case  fl3  is  a  triangle. 
My  (  t  .20),  we  see  t  hat 


/*;„(£)-(£  o’"  V'.(f)  (4.3i) 

where  y,  is  a  polynomial  in  £.  We  define,  for  0  <  s  <  (  —  1,  a  polynomial  by 

<;,„(£,»)  -  ({ '  ';.(())*;w(£4 »)'  h-32*) 

with 

^(o-  i).  (4.32b) 

t-n 

where  \*  1  and  for  s  0,  r;’m ( £ )  -  0.  We  now  define 

(4.33) 

•»-n 

Then,  using,  (4.29)  and  (4.32),  it  may  be  verified  that.  <ym  satisfies  (4.30a)  and  (4.30b). 
We  now  show  that.  (4.30c)  is  also  satisfied.  To  this  end,  we  first  show  that  r’m  has  a 
similar  decomposition  as  does  ft*m  in  (4.31).  Using  (4.32)  with  s  ==  0, 

o  _  /$»(o*»u  •  vY 
<,m  "  '  ((  iy 

so  t  hat  for  0  -  s'  f  I , 


nutw  !)■•■('  *  i  i) 

<hi'  U'  }  (£  - 1 )’ 

Using  ( 1.3 1 )  with  s  0,  this  yields 


'l •  .11 

•:  u. 

()ii 


i)  u  i )w  *  v;(o.  o-s 


(4.34) 


(4.33) 


f  I. 


Assume  mow  that  for  t  <  i  I, 


d'*j”  (£,~1)  -  tf  -  l)M-'-Vi(0  t<s<l-  1  (4.36) 

for  some  polynomial  <t>\.  Then,  by  (4.32),  (4.31),  (4.36), 

-- (t  i y-"mx,r(’m+r1)f 

lor  some  polynomial  (/>(£).  lienee,  for  i  <  s  <  t.  1, 


d'<\„ 

- ,m  (£, 
rltf  V 


1) 


u-i  y- 


*4 


<w 


di)’ 


(€.-') 


(4.37) 


for  some  polynomial  <f>\.  (4.35)-(4.37)  imply,  by  induction,  that  (4.36)  holds  for  all 
/  0,1 _ -  1 .  From  this,  we  see 


c(0-(e-  i)2,--v,(o 


(4.38) 


for  some  polynomial  <A«(£)- 

Now  by  the  definition  (4.32a)  of  we  have 


Since  on  S  —  T 


,  d*"*7  /M  '/ 

U  I 


C 

le 


vvr  srr 


Ikjm  lln'(.S)  <  ^  1)6^,  ^  (flj'mtt)  r/m(0)  il«"(r, 


I  dfi  -I! 


(4.39) 


I'siiiR  (4.31)  and  (4.38),  we  see  that 

hAUu,;-  r;"))(i1  " ,,>r 

Hence,  using  Lemma  4.2,  for  fs  f\  I  f2  )  f%. 


0, l,...,2 ( 


23 


I  </'■’ 


| ) r, )  L  '  7>"’)H,/"fr,)  -  C  (ll^ymll//'Mr.)  1  llr/mli///'.(r,))  H'10) 


We  now  show  that  for  i  —  0, ....  t  -  I ,  s  <  f  -  1,  r  <  f  . 

First,  we  see  that  using  (4.34),  (4.31),  Lemma  4.2,  and  (4.28),  (4.41)  is  satisfied  for 
i  -  0.  Next,  assume  that  (4.41)  holds  for  all  i  <n  -  1.  Then,  we  have,  using  (4.32) 


<9Yn 

lLm(c 

dn'  u’  ' 


so  that, 


)r(^*(0-^(0)U-l)<",+n 

Cl  n  1  dnC 

((  1)> "  £  olf 


<  e  I If  . . .  i  1  ...  'I  Iln-I.  ■■(!■,> 

(where  we  have  used  (4.31),  (4.37)  and  Lemma  4.2) 

<'  Cn-(k-f~r-  _ . 


by  (1.28)  and  our  hypothesis.  Hence,  by  induction,  (4.41)  holds  for  all 
i  0,1,...,/’  1 .  This  shows  that 


lr>m.lln''(r,)  Cp  '  f''  2 Ml ** II /f *(n) 


so  that,  using  (4.39),  (4.40),  (4.28)  and  (4.42)  with  Lemma  4.1,  we  see 


(4.42) 


ik;ji,r(„  '  c  y;  (P-i*  ,n|)(r'.-  !) 

t\  1  1  f.,<f 

'  Cp  {k  '’iMIn^n) 

Lsing  ( 1.33),  the  same  estimate  holds  for  ,  so  that  (4.30c)  is  proven. 

Hence  for  any  7  ft,  n  Um,  we  have  constructed  a.  polynomial  <,m  satisfying 

(1.30).  Defining  c;m  C we  see  that  ||cjrn  ll//'(n;)  also  satisfies  the  bound  in 


( 4.30c).  Moreover,  by  (4.30a.,b),  replacing  zj/1  by  zjd  c>m  on  tt}  achieves  (he  required 
(."  1  continuity  across  y  without  altering  the  jumps  in  zr  on  the  other  sides  of  d fl}. 
f?epeating  this  proress  for  eacli  y  in  the  triangulation,  we  obtain  a  up  *'  zr  £  C'{'  ’’(H) 
satisfying  (4.2).  The  essential  boundary  conditions  (4.1)  on  I'  may  be  imposed  on  up 
by  the  same  method.  This  completes  the  proof  of  the  theorem. 

□ 

Remark  4.1  Tlie  function  ur  constructed  by  us  belongs  to  0i>p+n  for  some  n  depend¬ 
ing  on  f.  Bv  suitably  changing  the  constant  in  (4.2),  we  obtain  ur  6  oSp  such  that 
(4 . 1  )-(4 .2)  still  hold . 

Remark  4.2  For  f  r  k  <  2 1.  *,  Theorem  4.1  still  holds  provided  we  assume  that 
u  c  <f>  instead  of  n  //^(H),  where  4>  is  defined  by  interpolation  using  the 

A-method, 

*  ---  (/r(n)  n  //'(n),  //'( n)),iCO  (4.43) 

Here,  r  >  2  f  |  and  (ffr  (Q),  Hf(U))ti_2  ~  Hk(Q).  The  proof  is  similar  to  that  of 
'I'heorem  4.2  in  |3)  and  is  omitted  here.  Generally,  however,  the  restriction  k.  >  2f.~  i 
is  not  a  severe  one,  particularly  in  the  light  of  results  in  the  next  section  where  corner 
singularities  are  treated. 

As  mentioned  in  the  introduction,  the  results  corresponding  to  Theorem  4.1  in  [9] , 
|IC>|  are  based  on  the  assumption  that  u  €  4’,  which  is  not  the  usual  result  predicted 
by  elliptic  regularity  theory. 

Theorem  4.1  and  Remark  4.2  lead  to  the  following  estimate  for  the  rate  of  con¬ 
vergence  of  the  p-version  of  the  finite  element  method. 

Theorem  4.2  Let.  u  f-  Hk (0)  n  Hf(Q),  k.  >  f.,  be  the.  solution  of  (2.2)-(2.f).  Assume 
further  that  for  f  <  k  <  2f  u  £  4>  defined  by  (f.fS).  Let  up  be  the  finite  element, 
solution  based  on  the  p-version  satisfying  (2.9).  Then 

llM  ”  Mrl|r/'(n)  S  ('P  ^  (4.44) 

where  (‘  is  a  constant  independent  of  p,  u  but  depending  on  the  partition,  of  fi. 

Proof :  The  proof  follows  from  'I'heorem  4.1  and  the  fart  that. 


5.  The  Approximation  of  Singular  Functions 


In  I  lie  previous  sect  ion,  we  analyzed  t  he  approximation  of  funct  ions  wliieh  were  known 
to  lie  in  n  k  >  2f  In  this  section,  we  analyze  functions  of  the  type 

(2.0),  which  have  a  singularity  at,  a  corner  of  the  domain. 

5.1.  An  Approximation  Result 

Let  Q  -  (  1,1)  x  (-  1,1)  as  before.  Let  f,  =  x,  \  1 ,  i  —  1 , 2  and  let  (r,0)  lie 

the  polar  coordinates  with  origin  (1,1);  r2  —  x2  I-  x\  ,  0  ~  arc  tan  (.xj/x, ).  For 

ft"  -  1 ,  0  p  <  1 ,  define 

~  {.t  Q  |  hr,  <  ,r2  <  ifi|} 

~  n  {.r  |  i2  I  .r2  <  p2} 

-  {.r  |  0  <  .r,  <  fl,  0  <  x2  <  «},  0  <  a  <  1  (5.1) 

Qh  -  {t  |  ,t ,  >  /),  x2  >  />}  H  Q,  0  <  b  <  1 

ft*  ■•■■■  .*vnQi 

Let  ft-0  I.  Fig.  5.1  shows  the  domains  under  consideration 


M,l)  (1,1) 


5S3 


^(.ri,.r2)  -  (x,  kx2)  {kx,  x2)'  -  rl  <f>{(6)  (5.2a) 

Obviously,  4>\{0)  is  an  analytic  function  in  0  and  £  is  a  polynomial  which  satisfies 

I)[k)f\  1.-0  0<\k\<?~\,  i  —  1,2  (5.2b) 

where  L\  and  I,2  are  the  lines  .rt  =  kx2  and  .f(  —  x.2f k  respectively. 

Let,  for  Re  n  >  f  -  I,  7  >  0  real, 

u{xt,x2)  =  Re  (r'Y|  log  r|1y(r)<^(0)}  (5.3) 

where  <h($),  \  (r)  are  sufficiently  smooth  real  (e.g.,  (7°”)  functions  and 

\(r)  -  I  for  0  5,  r  <  P 

,  2  p  1 

—  0  for  r,  0  <  n  <  - 

3  2 

is  a  function  defined  on  Q.  We  shall  assume  that  ?/  satisfies  (5.2l>)  on  1'uL  2  and  has 
support  in  UK„.  Then  we  see  that 


«(*l-*2)  £(.Ti,f2)un(x,,.f2) 


(5.4a) 


where 


"n(.r i,x2)  Re  {rn'  2' |  log  r|\\(r)'I'(0)}  (5.4b) 

with  ^(0)  a  smooth  (e.g.,  (7no)  function. 

The  main  approximation  result  we  prove  is 

Theorem  5.1  Let  11  be  given  by  (5.3).  Then  t here  exists  zp  f-  PP{Q)  such  that  for 
!)  \k\  f  1,  D^zr  -  0  on  the  lines  L,,i  —  1 , 2  and  for  k.„  >  k., 

II"  2(n<,a  'H)  (5.5) 

where  (!  is  a  constant  independent  of  p. 

We  will  require  a  series  of  lemmas  to  prove  Theorem  5.1. 

Let  '.i.'(r),0  r  -  00  be  a  f,"v'  function  satisfying 

w(r)  0  for  0  <  r  <  I 
1  for  2  <-  r  <  00 


UPw  (HVf 


WrW? 


For  anv  A  '•  0,  define 


Thru  \vr  decompose  un  by 


It  can  ho  readily  soon  that. 


u.-A(r)  w(  -■) 


t/n  ~  v  t  w 


v  —  oj  A  J/.,| 


w  -  (I  -  wA)u0 


v  —  0  for  0  <.  r  <  A 


w  0  for  r  >  2A. 


(5.8a) 


(5.8b) 


Lemma  5.1  (Uven  k  (/r,,/r2),  there  exists  a  constant  C(lc)  such  that  for  x  -■ 
{r\,.r2)  !<k„ 


\i)^v\  <  r;(/r) | log a|7(i  1  -w 
-  0 


on  IP 


where  d  Ro(o). 

Proof  :  For  o  real  the  lemma  follows  by  taking  or.  —  n  -  2?  T  2  in  Lemma  5.1  of  [3], 
The  result  for  o  complex  follows  easily. 


In  what  follows,  we  will  assume  that  v  satisfies  (5.9)  and  not  the  explicit  form 
(5.11)).  (5.8a). 

Let 


n(xi,.r2)  V^rt,;/.(.r,)/>(.r;) 


(5.10) 


i-11  y-n 


where  P,(rf)  P,{rf,ft,ft),  ft  -  I  are  Jacobi  polynomials  of  index  ft  which  will  be 
determined  later.  Then 


"o  f;.r; 


\<:A’  I  !)(>  I  I)  /  /  »’(^,,.r2)/’(.r,)/’J(.r2)(l  .T*)'Vr,f/.r2  ( 5.11) 


PI 


hiSvivIvS 


where  (  are  bounded  from  above  and  below  independently  of  i,j  but.  depending 
on  if  (see[ll|,  p.  811 ,  formula  7.391 . 1 ).  Define 


r  r 


i  -n  )-o 

(5.12) 

Mr  2) 

(5.13a) 

j-n 

(5.13b) 

w  i  1  h 

f».(r2)  f,(»  1  1 )/(  >’(r  l  •  rz)(l  ^)/,r,(T|)rfrl. 

(5.14) 

It  ran  lie  readily 

seen  t  hat 

1  n 

(5.15a) 

Let 

r  n 

(5.15b) 

t/v(rl-r2)  X^/,i(r2)/><(rl) 

1  -fl 

(5.10) 

l  ben 

»’  f’p  (»’  */’;■)  *  (V’p  »’p)  °V  * 

Pr 

(5.17) 

The  following 
and  5. 1  of  |3| . 

lemma  follows  immediately  l*y  taking  0 

0  2 f  t  2  in 

Lemma  5.3 

Lonmin  5.2  Let 

n  m  1  (!  2f  1  ]■  0.  Thru 

^  ’’I  ln*  A|7(1  1  ri)" 

*  '[  v*  \ 

(5.18a) 

,<1’"l’,{r2)  <'\  log  A|7  ,  m),; 

1  <ir7  1  (,  ,  ,)!  (’  '  T2) 

2f\  \ 

(5.18b) 

I,r1  iis  now  analyze  pv  (/>  given  in  ( .r> .  1  7 ) .  Wo  have 


Pp[t i,.r2)  '  XlM**)  /,!'l(-r2)!,’(r>) 

i-n 

so  I  Hal  for  k  (A-|,  fc2), 

1' ii..  i>!,'ii(x,)/’,i*,v1i 

t-fr,  "-r2 

Now,  for  /  -  0  integer,  wo  have  by  (5),  formula  A. 2. 15 

Pt|()(z,  /?,/?)  -  /- (2)3  I  *  4  1 )  . . .  (2/?  I  t  I  /)/>,  ,(.r,/?  I  /,/?  l  I) 
llonro,  for  0  ■  rn  4  k7  <  />  I  I , 


ii  f  ii 


U?(K  ^)(x,))  ° 


<  I  j  ;■  i  I 


r  v,„,  ,  y  1  m  1  ‘-4' 

i';"  r;  ,y,  U  ’»  h>'i 


■  (I  (*■«») 

I  sing  (5.1  1),  wo  soo  tlial  tho  support  of  b,(r7)  lios  in  lt  |  I  I  Asinfln,0|  wlioro 
tanff,,  1  .  llonro,  from  (5.181))  and  (5. If)),  for  o  rn  I  f  2f  I  J  •  0, 


<  ‘  r  f"  i  , '  *■-’/»  \  2 

7  .  c 

n/m  y  i,A  *m«.,  \  ff.r  *  1 


<:\ log  Ap 


r  /-ii 

7. 


(  I  A  ini  I1,,  f 


(I  i  .r2) 


2-.  m  i  ? a  «r  *•»  i 


•”  !rf.r, 


4  |  log  Ap 1  2.»  ™u/i  if  *  ; 


,,2m  2fri  I  S 


(5.20a) 


for  (ho  case  'l,  provided  2<\  m  I  2  ft  If5  k 2  1  2  '  0.  For  U10  raso  k j  0,  wo 

use  (r».jSa)  instead  of  (r>.|8h)  to  hound  '*  and  obtain 


I  |  log  A  |21  2-»  ml  2/1  M  k ,  I  '• 

11  '  J 


(5.20b) 


provided  2 A  rn  *  ‘2/^  •I/’  /r2  I  '  0. 

Similarly,  we  estimate  the  term  I)^ar  with  0  '■  rn  \  k2  '  p  t  1 ,  rn  "•  k\.  We  have 

*'>n 

rrr(-rt ,  -r2)  l\(r2)P,(rt) 


/?.  /  '  ( 7  ( ’  ( /2|ir,rrr  )2(  I  rjrS/r.Hl  .r2)" 1  ^r/.r2 


.rH)',,*Jr/.r2 


f  /^/>,(.r2)\2  (»  I  m)! 
‘■|l  '  i  ,  mi  V  ^r2J  /  »(* 


i  /■"  /  /)”"*? i 


■  •,!./  ,  /  , 


<)tT<)t27 


r2)^  1  mr/.r |  ( 1  .rjr^rj  (5.21) 


Since  tho  support  of  r  lies  in  //„n  .S'^  ,  we  can  use  Lemma  5.1  anrl  obtain  with 

/>  1  i  1  (I  I  X2).  I  l  K„(l  ,r2)| 


( 'j  log  A| 


f  /" 

11  **  I  I  A  im  K„  I- 


2'  m  ^>(|  T2)/,,n’(l  .r2)/,,'r3J.rl(V.r2 


1 \  1  re 


r'S  log  A 


vr  [ 

M  J 


2(m  •'ha  «m  0, 


(I  t  ,r2) 


2->  m  f  2/1  4'  »•;  t  1 


(  !  log  A  I  1  2,y  ,  J,,  Ir2  I  2 


.2  (m  k  i  | 


(5.22) 


provided  that  2  A  rn  l  2/f  4f  k7  I  2  '  0.  lienee,  we  obtain  the  following  lemma. 

Lomm.')  5.3  l,rl  pv  ond  ny  hr  ns  defined  in  (5.1 7).  1  hen  for  0  <  w  I  f  p  I  I  nnd 
0  A- !  f 


/ ' 1  / "(o"VrP(i  i  iji'"’*i*i  • 


where  k\  maxjl,^}  and  f  ~  I  if  /rt  —  0,  f  0  otherwise,  provided  that  ir  rn 

C  2f  !  ]  ■  0  and  2n  m  I  2/,  M  I  |  ^  0 


/  /  (/><MrT,,)2(  I  Ti ),?  1  *’  ( I  •r2)/,  *  k‘d.r\dx 2 


ciiogAr 

n2(»i  *•  i ) 


^2,<  r>,  12/1  if  *2  I  2 


(5.24) 

provided  that  2 1\  n>  I  2/1  \(  t  2  0.  The  constant  C  is  independent  of  k.A.p  hut 

drp'iids  on  (\.d,'),in,f. 

Let  Q 2A  hr  as  in  (.r).l)  and  define 

it*  nK  r  i  (y2A 

\lso.  for  /(•/■), .r2),  (.r,,.r2)  <  Q  and  A  define 

/a (  r i  >  -r2 )  /(r,  2A,  .r 2  2A).  (.r:  ,r2)  r:  y2A 

0  (.t,,t2)  r  y  y2A  1  j 

Loimnn  5.4  Le/  £(r|.;r2)  he  given  hy  (5.2a)  and  let  0  <  A  <  Then,  on  /?An, 

i/^(Ma(T,,r2)|  <.  <7(1  x?)'((l  .r2)'3  (5.20) 

for  uny  t  j,/2  •  0  satisfying  I,  I  f;  2f  |  /r  | . 


(5.25) 


Proof :  The  proof  is  essentially  thr  same  as  that  of  Lemma  5.5  in  [3|  ami  is  omitted 
here. 


Lemma  5.5  Let  »•  satisfy  (5.9)  and  r,,  he  given  hy  (5.12).  Then  for  A  -  p 

IIsa('’  »v)IIh'(/?ai)  •  < '| Ior/»|>  7h  n  11 
where  <7  is  a  roust  ant  independent  of  p. 


(5.27) 


Proof :  W r  first,  estimate  || II // '( I  where  pr  (and  or)  are  as  in  (5.17).  To  this 
end,  let,  us  estimate  f)t  -  |l(^*r*^A)( )  with  |r|  I  |.s|  (.  Using  (5.25), 

we  have  for  any  U,f2  0  with  f.\  I  /2  2 f  |r| 
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\\r  may  similarly  show  that  for  |r|  (  |.s|  f. 
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I  h  is  proves  t  hr  lemma . 


fur  the  ease  that  d  '  I,  d  not  an  integer,  let  k  -  [d|  he  the  largest  integer  less 
than  d  l  or  7  an  integer.  0  -  7  k,  wr  denote  by  t>m  the  7th  derivative  of  v  along 
the  direction  n ,  where  i>  is  the  unit,  vector  along  t  lie  line  T|  t2.  Then  will 

satisfy  (",.«>)  in  Lemma  .">.1  with  d  replaced  by  d  7  '  0.  Hence,  using  Lemma  .r>.5, 
we  get 

i|£A(»‘w  r;.'!,!l»',wA  ,  r |l«gp|>  2'"  7  (r»  28) 

Let  1,0  defined  by  (.'».(>)  and  u’1}  be  its  translation  given  by  (.r».2.r») ,  Then  (see  (  V7)) 


"r>A  ttnAwA  *  ’'"aU  1 

»’A  I  '"A 

l'»in<  <■  »  1  Uf{U,„),  Hien  ?/A  ‘  // ^ ( /»'*,, )  and  hence 

£a«*a  «a(i  //W 


(S-2») 


iV»V.i'iW 


I.eiimia  r»,0  l.fl  A  p  A  2\  2A.  Thru  for  k  Ini,  n  f 


k 

\  'i 


A' 


<-\WrVr 


2(,. 


i|£.\ "’a il//,(/^ii)  ■  < ’I  l 2("  " 11 

ulirrr  ('  is  iiuirpruiirnl  of  ]>,  A  . 
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In  tlir  above  inerpia lit v,  wo  used  the  obvious  farl  that  o  k  f 
tin1  other  terms  in  (o. .'!())  analogously. 
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wIhto  wo  havo  used  I lio  fact  dial  d  ■  f  I.  'I'lio  other  terms  in  (5.31)  ran  ho  treated 
similarly.  This  completes  the  proof  of  Lemma  5.6. 
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We  now  prove  our  main  result. 

/ 'roof  of  'theorem  5,1 

hot  r,  r,  2A,  >  1,2  and  let  SK  he  the  translation  of  SK  obtained  by  this 

transformation.  Let 
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where  we  have  used  (5.28).  (5.30)  ami  (5.3|).  We  now  translate  hark  to  SK  and 
s'lil  al>U  adjust  the  constant  ('  in  (5.5)  to  gel  the  theorem. 
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Remark  5.1  We  have  proven  a  slightly  stronger  result  than  Theorem  .r».|.  It,  is 
siiUi<  dent  to  assume  that  v  and  t r  defined  hy  (5.Ra),  (5.8b)  sat.i  fy  (5.9)  and  (5.32) 
respect  ively. 

Uoiiiark  5.2  lor  the  rasp  that  lira  is  an  integer  and  7  ~  0,  ti  will  he  arbitrarily 
smooth.  Ilenre.  the  above  result  is  too  pessimistic  and  the  results  from  Section  4  will 
apply. 

Remark  5. ."4  From  the  proof,  it  may  be  seen  that  the  internal  angle  w,  between  7, 
and  ,  1  could  equal  2tt,  i.e.,  we  may  also  consider  the  slit,  domain. 

5.2.  Approximation  over  the  domain  U 

We  now  return  to  the  problem  of  approximation  of  the  functions  j/.^  given  by  (2.6). 

bet  the  vertex  .4,  of  U  be  at  t he  origin  0.  bet  the  part  of  0  containing  the  elements 
with  vertices  at  0  be  as  shown  as  in  Fig.  5.2.  We  assume  tht  we  have  only  triangular 
elements.  The  rase  when  elements  are  parallelograms  does  not  change  the  argument. 


Tm  V 


Figure  5.2 

bet  \)  '  i"’  ,7,,  I’  1  'j"  ?  H,  H, ,  | .  bet  the  line  0/1,  have  t  he  coordinate  # 

I .....  m  1  1 .  Denote  />(,  {.r  |  r]  I  x\  ■  p }  and  assume  that  (  U,  0  -  f 

We  then  obtain  the  following  theorem. 


Theorem  5.2  Let  ?/  hr  the.  function  given  by  (5.3)  with  p  pnp,  (i  sufficiently 
small.  then  there,  exists  zr  f-  satisfying  zv  f  Pr(7\ ),  t  ~  1 ,  /7^2p  0 

on  ()li\,()J fm  )  |  rr r? r/  1'  for  0  <  |r|  <.  f  I  om/ 

II w  *rll/f'(n)  r^’|loR,,|>  S<B"*  'l|>  (S. 33) 

where  < .'  is  independent  of  p. 

The  proof. of  the  above  theorem  is  very  similar  t,o  that,  of  [ 3 ] ,  Theorem  5.2  and 
only  a  brief  outline  is  given  here.  Essentially,  we  first,  consider  the  case  for  which 
lfr]<\'(0})  0  for  ;  -  I  1,  0  <  jr|  <  f.  ■-  1.  We  may  then  map  S  — 

{(r.0)\0j  0  -  0; , ,  }  onto  IiK  by  a  linear  mapping  T  and  consider  the  image  rt  of  i/ 

on  Tj  T(Tj).  Let  t]}  be  a  polynomial  function  of  degree  <  f.  satisfying  D*r'r;;  -  0 
for  0  <  r  <  f  l  on  T [B^B,^  Then,  after  suitably  extending  u  outside  RK , 
the  function  u j 1]}  satisfies  the  conditions  mentioned  in  Remark  5.1  to  Theorem  5.1. 
Hence,  we  may  approximate  u,jp,  by  a  function  z"v  satisfying  (5.33)  on  7)  and  hence 
Zj.1],  z‘f  satisfies  (5.33)  too,  proving  the  result  for  this  case. 

For  the  case  when  /7*r,?77  —  0  for  j  f  /e,  0  <  r  <  t  1,  the  triangles  T;„  .j,  7’;„ 
are  mapped  together  into  fiK  and  the  argument  repeated.  The  details  mav  be  found 
in  [ 3 1 . 

Remark  5.4  The  function  we  constructed  was  in  Pru(7])-  Hy  suitably  changing 
the  constant  in  (5.33),  we  may  obtain  a  function  in  Pn( 7\). 

F?oiiinrk  5.5  (5.33)  obviously  yields  the  estimate 
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0.  The  Hate  of  Convergence  of  the  p- version  of  the  Finite 
Element  Method 

Wp  now  summarize  our  results  from  Sections  4  and  5  and  briefly  remark  on  some 
generalizat  ions. 

The  following  theorem  follows  immediately  from  Theorem  4.2  and  Theorem  5. 2. 

Theorem  G.l  I,et  u  be  the  solution  of  problem  (2.2)-(2.f) .  Assume  that  u  ran  be. 
written  in  the.  form  (2.5),  (2.6)  and  in  addition,  that  for  f.  <  k  <  2(  u  f:  <J> 

defined  bp  (f.jX).  Let  ur  be  the  finite,  element  solution  as  described  in  Section  2.5 
with  triangular  or  parallelogram  elements.  Then 

II"  -  "t II //'(d)  <  Cp  “\\ogp\t'R  (6.1) 

where,  letting  n‘  ~~  minReo^, 


//  min(A:  <!, 2(R.eo,  t  -I  i))  ■-  min(/r  -  i,  2(a*  -  (. -\  1))  (6.2) 

i >  ~  max{7^  :  Re  rv^  -  o*}  if  p.  —  2(o'  -  t- 1-  1) 

-  0  otherwise  (6.3) 

ft  ■■  il "  i  II //‘{(i)  i  K'!'1! 

') 

Remark  0.1  'I’heorem  6.1  has  been  slated  only  for  the  model  problem  (2.2)-(2.d). 
It  is  obvious,  however,  that  the  theorem  holds  Tor  any  elliptic  problem  or  order  2 f 
if  the  solution  has  the  form  (2.5),  (2.6)  or  when  (2.6)  is  different,  but  has  the  same 
rhararter  concerning  the  growth  of  its  derivatives.  Moreover,  as  mentioned  in  Section 
2.1,  more  general  boundary  conditions  may  be  also  be  treated. 

Remark  0.2  We  assumed  that  polynomials  of  the  same  degree  are  used  over  each 

element.  Our  results  and  proofs  may  hr  modified  in  an  obvious  wav  when  different 

degrees  are  used  over  different  elements. 
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